Dynamics of fluctuations in an optical analogue of the Laval nozzle 



O 

(N 

(N 



C/3 

o 

• I— I 

> 

Oh' 

q : 

o ■ 
•I— I ' 

Oh 



(N 
> 

m 
m 



On 
O 



X 



Itzhak Fouxon/ Oleg Farberovich/ Shimshon Bar- Ad/ and Victor Fleurov-'^ 

^Raymond and Beverly Sackler School of Physics and Astronomy, Tel- Aviv University, Tel- Aviv 69978, Israel 

(Dated: September 27, 2010) 

Using the analogy between the description of coherent light propagation in a medium with Kerr 
nonlinearity by means of nonlinear Schrodinger equation and that of a dissipationless liquid, we 
propose an optical analogue of the Laval nozzle. The optical Laval nozzle will allow one to form a 
transonic flow, in which one can observe and study very unusual dynamics of classical and quantum 
fluctuations, including an analogue of the Hawking radiation of real black holes. Theoretical analysis 
of these dynamics is supported by numerical calculations, and estimates for a possible experimental 
realization are presented. 
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Black hole radiation is one of the most impressive phe- 
nomena at the intersection of general relativity and quan- 
tum field theory. Accounting for the quantum nature of 
the physical vacuum led to the prediction that a black 
hole, defined classically as an object that even light can- 
not escape, in fact can be characterized by a temperj,- 
ture, entropy[H and moreover emits thermal radiation[2]. 
Since direct experimentation with black holes is hardly 
possible, it was suggested to consider analogous phenom- 
ena in condensed matter physics where the "high-energy" 
(short- wavelength) physics is known The suggestion 
was based on the observation that the derivation of the 
Hawking radiation uses only the linear wave equation 
in curved space-time, and not the Einstein equations. 
The same conditions for wave propagation arise when 
considering sound propagation in a fluid when the back- 
ground flow is non-trivial This similarity was shown 
to be sufficient for applying the original considerations 
by Hawking, predicting the existence of thermal radia- 
tion of quantum origin from the fluid counterpart of the 
horizon, which can be called Mach horizon (where the 
fluid velocity equals the sound velocity, i.e. Mach number 
M = 1). These ideas were further developed forpossible 
applications in BEC fluids and other systems H-Q. As 
for experimental realizations, a white hole horizon was 
observed in optical fibers |9j, where the probe light was 
back-reflected from a moving soliton, and a black hole 
horizon was observed in a BEC system[l3. 

We propose an experimental setup capable of creating 
a Mach horizon in an optical medium with Kerr nonlin- 
earity, which may be called optical Laval nozzle. The 
propagation of weakly nonlinear coherent optical pulses 
is described by the non-linear Schrodinger (NLS) equa- 
tion 

r>A 1 

^%^ = -77^^^A + U{x,v)A + \\A\''A. (1) 

It assumes the paraxial approximation, according to 
which the electric field of the light wave is written as 
E — A{x, y, t; z)e~'^°^, where A{x, y, t; z) is weakly z de- 
pendent complex amplitude of the light propagating in 



the z direction. The time coordinate t is converted into 



the r coordinate [Ul which describes the shape of the 
light pulse in the moving coordinate system. The radius 
vector differential is now dr = {dx,dy,dT) (anomalous 
dispersion) or dr = {dx,dy,—idT) (normal dispersion). 
U{x,y) is the equivalent external potential created by 
spacial variation of the refraction index in the medium 
and does not depend on t. 

Now we briefly sketch the fundamental result Q as ap- 
plied to coherent light propagation in a medium with 
Kerr nonlinearity. The Madelung transformation 12 
A = /e^"^ (see also[iMi) maps the NLS eq.(II|) on two 
equations 



-f V • [pv] = 0, 



/3o 



(2) 

V {Vq^ + U + \p) (3) 



for an equivalent fluid with density p = and velocity 
f3oV — —V4>. For permanently shined light the coordi- 
nate T is redundant. Vqu — ~2^^j^ is a "quantum 
potential" corresponding to h = 1. The light wave vector 
j3o is the " particle mass" and the velocity v is dimension- 
less. 

Linearizing Eqs. ([2]) and ([3|) with respect to the small 
fluctuations p—po = poip and 4' — 4>{) — </? around a steady 
solution po{x,y) and (po{x,y), we arrive at 



(4) 



where the metric g'^'^ with the determinant g is given by 
the interval 

da^ = g^^dx^dx'^ = 




dr — 2dzvQ ■ dr 



/3o 



dz' 



(5) 



with g^'i/ — Sj^. The ordinary "sound waves" in the 
effective fluid arise as solutions of eq. (|4]) around the equi- 
librium solution pq = const, (f>o = const and vq = 0, that 
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exists at [/ = 0, so that Xpo//3o = is the squared sound Equation ([2]) for a steady flow becomes 



velocity. The nonlinearity coefficient is positive, A > 0, 
otherwise the "sound velocity" becomes imaginary and 
various instabilities, such as collapsing solitonsfll] arise. 
Eq. (HI may include corrections due to the quantum po- 
tential, which may be of importance for short waves. 

The background flow may arise in the Laval nozzle 17 , 
fisj l. i.e. a vessel with variable cross-section ^(a:), whose 
application to some condensed-matter analogues was dis- 
cussed in[l^ E^l- The flow velocity increases with de- 
creasing S (cc) , and reaches the sound velocity in the nar- 
rowest part of the vessel, called the throat. Further ac- 
celeration is reached by increasing S{x). Here we ana- 
lyze the transonic flow for the optical realization of the 
Laval nozzle shown in Fig. ??. The laser beam, nearly 




FIG. 1: A cartoon of a laser beam propagating inside the 
optical Laval nozzle. The beam, with a small initial angle 
relative to the z axis, bends more towards the x direction as 
it propagates along the z direction ("time"). The bend in the 
cartoon is strongly exaggerated. A fluctuation is schemati- 
cally shown near the Mach horizon as it is cut in two parts 
which propagate in opposite directions. 

parallel to the z axis, is tilted in the x direction by an 
angle arctan(fcj;//3o), i.e. the initial phase in the elec- 
tric field amplitude A{x,y) is iy9 = k^x. k^/Po plays the 
part of "flow velocity" in the x direction. Acceleration 
of the flow within the plate of variable width (y - axis) 
results in bending of the beam {k^ increases) as shown 
in Fig. ??. It leads to small corrections to /3o, quadratic 
in kx/Po ^ 1. We always assume that the wave length 
1/Po is smaller than all the other relevant scales. 

Now we consider the flow bounded by the walls of a 
hyperbolic shape 



X 
62" 



1 



(6) 



with sliding boundary conditions, and the elliptic coor- 
dinates X — csinhui sinu2, y = — c cosh mi cos M2, are 
introduced, with c = ^/aF+lP. From the symmetry con- 
siderations the flow is along the x axis for y = 0, i.e. 
U2 = 7r/2. 



d d 
— VA(pwi) + — VA(pw2) = 

OUx OU2 

where A — sinh^ ui -t- sin^ U2 and vi and V2 are the ve- 
locity components in the elliptic coordinates. 
This equation is satisfied if 



pvi 



PV2 



Fi{u2) + fui 



VA 

F2{ui) - fu2 



(7) 



where -^1(1*2) and F2{ui) are functions only of U2 and ui 
respectively, and / is a constant. Since the transversal 
velocity V2 is zero at the throat at ui = we get that 

/ = o. 

We assume the streamlines in the vicinity of the central 
line, i.e. small U2 — U2 , to be close to the lines of 

constant U2. Then calculating the covariant derivatives 
in the curvilinear elliptic coordinates and neglecting in 
the result the second order corrections in ■U2, Eq. ^ 
reads (cf., [11]) 



djpvi) 
dui 



Xp 



dvi 

dui ' 



(8) 



The flow reaches the sound velocity at ui = 0. The 
second Eq. ([7]) in the elliptic coordinates results in 
p{ui)vi{ui) — d — \dul with d = Fi{tt/2) = p s, and 
finally, using Eq. ®, we get vi — s{l + p = 

p(l — where overlined quantities refer to the throat. 
The coordinate dependence of the density immediately 
yields the coordinate dependence of the sound velocity 
\p{ui) = /3os^(ui). Correspondingly the values of these 
quantities at the throat are connected as Xp = Pos^ . 

We may now represent Eq. Q in the narrow transonic 
region in the dimensionless form 



2d^d(;(p + aS,d'^(p + aid(;(p + ad^(p — 0, (9) 



where C = —z, f = ui « — . The coefficients ai = 
c c 

l/VS, and a — -s/S are obtained above for the hyperbolic 
Laval nozzle. However, the ratio a/ai — 3 holds for any 
differentiable shape of this type. We look for a solution 
of Eq. (??) in the form 93(^,0 = exp{i;/C}^(^), where 
the functions ^(^) are linear combinations of the two sets 
of eigen functions 



and </J2(0 



(10) 



In a background flow characterized by a space indepen- 
dent velocity, eq. ^ has the form of a one-dimensional 
Klein-Gordon wave equation for the field ip ina flat space. 
Then one obtains two plane waves propagating to the left 
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(against the flow) and to the right (with the flow). In the 
considered here transonic case, when the background flow 
velocity and speed of sound vary with the coordinate, the 
left moving plane wave transforms into the eigen mode 
which is a function with a branching point near 
^ = 0, whereas the right moving plane wave is character- 
ized by the nonlinear spectral relation 



(11) 



The complex wave vector for the real "frequencies" re- 
flects the fact that the coordinate ^ of these eigen 
modes corresponds to motion along diverging hyperbolic 
streamlines. However, the frequency v is real which in- 
dicates linear stability of the flow with respect to weak 
fluctuations. 

The corresponding eigen modes 



V'i(0 



s 1 + 

/ 



^i(c), 



cp 



(12) 



for the density fluctuations can be obtained directly from 
the linearization of the Euler equation ([3]). 

Crossing the Mach horizon at ^ = 0, the solution 
^i(^) acquires the factor exp{^}. The same factor ap- 
pears when considering quantum fluctuations. Following 
the logic of j3| (also discussed in detail in [2l'| — eqs. 
(2.88)-(2.90)) this results in radiation to the left of the 
Mach horizon (outside the "black hole") with the fre- 
quency distribution N{v) = l/(exp{^} — 1). This dis- 
tribution has the same form as a Planck distribution for 



black body radiation, with the parameter Th 



physical dimensions. The latter plays a role similar to 
that of the Hawking temperature in black hole radiation. 
However the quantity Th is not a real physical temper- 
ature. The part of time is played here by the spatial 
coordinate z, i.e. the "frequency" v is the wave vector 
of a fluctuation in the z direction. This "temperature" 
is measured in units of momentum rather than energy, 
and the corresponding "thermal fluctuations" are waves 
along the z direction with a typical wavelength of order 
Xh ^ h/Tn ^ 2ttc/{V3s). 

A "straddled" fluctuation 



exp {i [i^o + \/3rV0(e)] (\/3 In + C) } 

exp 



V3 



(13) 



{0{£,) is the step function) with a Gaussian spectral 
density g{iy) = ^^P *''^2r^^ } S'i'ound a pos- 

itive frequency vq of a width T is composed exclu- 
sively of the "left moving" normal modes ^i(0: ^iH" 
gular at ^ = 0. Its "time" evolution {i.e. ( depen- 
dence) is determined by the envelope function G(^, C) = 

exp (VSln 1^1 - C)^} • The "frequencies" of these 

oscillations on both sides of the horizon strongly vary 
with ^. This is reminiscent of "time slowing" near the 
horizon of a real black hole. 

Fig. [5] shows how the maxima of the fluctuation ampli- 
tude propagate on both sides of the Mach horizon, with 
zero amplitude on the horizon. The part residing to the 
right side of the horizon (the supersonic region, ^ > 0) 

should be multiplied by the factor exp | ^ + ^^-^ | > 
1 (not shown in Fig. [2]), with the Hawking temperature 
TTVo/VS = hh'o/2TH as a parameter. Therefore the ratio 
of amplitudes of a classical fluctuation on the left and 
the right of the horizon directly measures the Hawking 



temperature. 

The fluctuation moves against the transonic flow with 
the velocity of sound. Its left part escapes to the left 
against the subsonic background flow, and moves with a 
very small velocity s; — u/ ~ s — to the left. The right 
part of the fluctuation is "flushed" away by the super- 
sonic flow (see Fig. ??), and moves to the right with a 
small velocity Vr — Sr ^ Vr — s. Here vi and s are the sub- 
sonic flow velocity and local speed of sound of the flow in 
the vicinity and to the left of the Mach horizon, whereas 
Vr and Sr are the corresponding quantities in the super- 
sonic region to the right of the Mach horizon (Note that 
for the symmetrical potential and spectral content of the 
fluctuation studied here the propagation on the two sides 
of the horizon is approximately symmetrical, as shown in 
Fig. [2]). Using the eigen modes Q for the density fluctu- 
ations we get a similar picture. In the case of quantum 
fluctuations the part moving to the left is analogous to 
the Hawking radiation observable outside the black hole. 
The part moving to the right is not observable in a real 
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FIG. 2: Propagation of the straddled wave packet to the left 
and right of the Mach horizon, = 0, (2 = 0.5, 1^3 = 1, 
C4 = 2 with r = 2. 



black hole. However in our case the supersonic region is 
accessible for observation. 

The normal modes {^) in (1101) are plane waves propa- 
gating to the right in the direction of the flow with the ve- 



This velocity is twice the sound 



locity V2 = s(^-^ 
velocity, i. e. the sound velocity on the background of the 
supersonic flow, for large enough " frequencies" ly ^ 1. At 
small frequencies ^ <C 1 the propagation velocity is closer 
to thrice the sound velocity, i.e. twice the sound velocity 
relative to the background flow. If we use the analogy 
with the real black holes this result would correspond to 
superluminal motion. However there is no paradox here 
since the wavelength is larger than the throat width in 
the Laval nozzle, 1/fc ^ c (in physical units). Hence we 
deal with the narrow fleld effect, analogous to superlumi- 
nal propagation in near fleld optics and tunneling [22 - 24 1 , 
which does not violate causality. 

The above results are illustrated by numerical simu- 
lations of the NLS eq. (HI in 1 + 2 dimensions using a 
finite-difference approach based on the split-step Crank- 
Nicholson method. [25] The potential U{x,y) is zero out- 
side the hyperbola (|6]) and a negative constant Uq inside. 
A packet with a Gaussian density distribution is given a 
Galilean boost in the fashion similar to Ref.[26| with a 
velocity close to that of sound for the density in the maxi- 
mum of the Gaussian. This configuration differs from the 
discussed in the analytical part of the paper and is chosen 
for its relative simplicity and due to the fact its exper- 
imental realization seems to be rather straightforward. 
Although it differs in many aspect such as strongly non- 
homogeneous distribution of the density and, hence, of 
the sound velocity, which is also affected be the quantum 
potential contribution, this configuration may be very 
convenient and useful for studying propagation and clas- 
sical fluctuations in the optical Laval nozzle. 

The calculation is carried out for Eq. ([T|) with /3o = 5 



and A = 10, whereas the amplitude function A is normal- 
ized to unity. The depth of the potential in the nozzle 
Uo — —200. It has a hyperbolic shape with a = 0.3 and 
b ~ 7.5. Meaning that its half-width is 0.3 at the throat 
at X — and 0.5 at x ~ ±10. The Gaussian centered ini- 
tially ata; = — 1, j/ = Ois given a boost with the phase 
if — 2.5x corresponding to the unit velocity. Figs. [3] 
show the packet at four successive stages of its propaga- 
tion. A straddled fluctuation gradually develops with a 
characteristic dip close to the throat of the nozzle, which 
becomes deeper and nearly cuts the original packet in 
two. The fluctuation in this case becomes strong and in 
the end consumes a major part of the packet. 

Since the numerical code deals with the NLS eq. ([T]), 
it solves a nonlinear problem strongly deviating from the 
steady transonic flow assumed in the above analysis, and 
fully accounts for the quantum potential. The density 
is strongly nonhomogeneous which leads to a nonhomo- 
geneous distribution of the sound velocity. The latter 
should be also corrected by the contribution of quan- 
tum potential which cannot be considered small any- 
more. The Gaussian packet evolves into a two hump 
structure permanently changing its shape with the prop- 
agation length z ("time"), as is clearly demonstrated in 
Figs. El 

There are such factors as quantum potential (disper- 
sion) and defocusing nonlinearity, which may influence 
the dynamics of the packet and in principle also cause 
its splitting in two parts. Control calculations of the 
packet moving in free space, or in a potential well with 
parallel walls do not show the type of behavior shown 
in Fig. 131 Decreasing the boost velocity also weakens 
the effect. Otherwise the general pattern presented in 
these figures is rather robust and holds when varying the 
parameters as long as the packet is boosted in a hyper- 
bolically shaped potential well.F 

The above ideas can be tested experimentally by 
launching a continuous wave laser beam through an ap- 
propriately shaped nozzle, with refiective walls, filled 
with a Kerr-like defocusing nonlinear material. Such a 
one-dimensional nozzle can be constructed from two con- 
vex cylindrical mirrors, put back to back at an adjustable 
separation with the space between them sealed and filled 
with a nonlinear liquid such as iodine-doped ethanol 27 1 . 
The nonlinear coefficient A of the solution can be ex- 
pressed, in terms of the nonlinearly-induced refractive 
index change (5n, as Sn0o/no. where no is the linear re- 
fractive index of the material[27|. The corresponding di- 
mensionless sound velocity is s = ^ 5n/no. Assuming a 
nonlinear index change of 1%, an incoming laser beam 
would emulate supersonic fiow when k^/Pa > 0.1, i.e. 
when the input angle of the beam is larger than « 6 
degrees. The required input angle for observing tran- 
sonic acceleration through the nozzle is therefore tech- 
nically feasible. The acceleration can be controlled by 
mechanically varying the throat width, and monitored 
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z =^ 





FIG. 3: (Color online) Four figures (a) to (b) show snapshots at different stages of the Gaussian packet motion through a 
hyperbolic Laval nozzle and development a straddled fluctuation with a characteristic dip close to the throat of the nozzle. 
The arrow in figure (a) to (d) shows the center of the packet if it were moving in a free space, which corresponds to the shift 
Ax = 1, 2, 3 and 4, respectively, from the starting point x — —1. The propagation distance z ("time") is indicated at each 
figure. 



by observing the deflection of the laser beam as it tra- 
verses the nozzle. The intensity and phase profile of the 
beam coming out of the experimental setup can be stud- 
ied using standard detection and interferometric tech- 
niques. These vi^ould allovif measurement of the transverse 
acceleration of the beam, as well as detection and char- 
acterization of classical fluctuations that propagate away 
from the throat region, such classical fluctuations may 
be introduced by locally perturbing the optical fluid in a 
deterministic fashion, e.g. by injecting a second narrow 
laser beam into the nozzle throat [l5j|. 

Constructing a Laval nozzle of hyperbolic shape might 
not be easy, but we believe that most of the effects 
discussed above will also be observed in other types of 
convergent-divergent nozzles (in fact many aeronautical 
applications of the Laval nozzle are not hyperbolically 
shaped). Otherwise the proposed experimental setup is 
quite simple and straightforward to implement in a table- 
top experiment, in particular as compared to the previ- 
ously proposed rotating black hole configuration 15| and 



to water tank experiments [28| . The (classical) fluctua- 
tions will be observed directly, superimposed on an other- 
wise smooth transverse beam proflle, and in principle will 
not require as high a dynamic range as previous optical 
experiments [9|. The experiment is also more straightfor- 
ward than other proposed schemes, e.g. those based on 
SQUID array transmission lines Q and atomic BEC[^[3]. 

In summary, we demonstrate the possibility of creating 
an optical analogue of the Laval nozzle. The interesting 
point will be to study the dynamics of straddled fluctua- 
tions which may be either quantum or classical and even 
artificially created. The equivalent of the Hawking tem- 
perature enters as an important parameter characterizing 
all types of fiuctuations. This temperature is measured 
in units of momentum, rather than energy, which cor- 
responds to a wavelength exceeding the width of nozzle 
throat by about two orders of magnitude, and in princi- 
ple accessible for experimental measurements. Numerical 
simulations support the theoretical findings but also in- 
dicate that the phenomenon is much more complicated 
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and interesting. 
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